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ABSTRACT 

Recent observations have demonstrated a significant growth in the integrated stellar mass of the red sequence 
since z = 1, dominated by a steadily increasing number of galaxies with stellar masses < 10 u Mq. In this 
paper, we use the COMBO-17 photometric redshift survey in conjunction with deep Spitzer 24/im data to explore 
the relationship between star formation and the growth of stellar mass. We calculate 'star formation rate functions' 
in four different redshift slices, splitting also into contributions from the red sequence and blue cloud for the first 
time. We find that the growth of stellar mass since z = 1 is consistent with the integrated star formation rate. 
Yet, most of the stars formed are in blue cloud galaxies. If the stellar mass already in, and formed in, z < 1 
blue cloud galaxies were to stay in the blue cloud the total stellar mass in blue galaxies would be dramatically 
overproduced. We explore the expected evolution of stellar mass functions, finding that in this picture the number 
of massive > 3 x 1O 1O M blue galaxies would also be overproduced; i.e., most of the new stars formed in blue 
cloud galaxies are in the massive galaxies. We explore a simple truncation scenario in which these 'extra' blue 
galaxies have their star formation suppressed by an unspecified mechanism or mechanisms; simple cessation of 
star formation in these extra blue galaxies is approximately sufficient to build up the red sequence at < 10 1 Mq . 

Subject headings: galaxies: general — galaxies: evolution — galaxies: stellar content — infrared: galaxies 



,4 



o 
o 

u 

Oh: 

< 

m 

6 



> 

o 

m 

o 
o 



X 



1. introduction 

Understanding where and when the stars in galaxies form 
is one of the key goals of observational extragalactic astron- 
omy. In the last decade, magnificent progress has been made 
towards this goal, with the construction of large galaxy red- 
shift surve ys allowing the evo l ution of cosmic star formation 
rate (e.g.. iMadau et al] 119961: iLillv et all 119961: ISteidel et all 
1999 ; iFlores et all 1 19991: iHaarsma et al.l l2000b iHopkinsI 120041 



Le Floc' h et all 20051) and the build-up of stellar mass (e.g., 



Brinchmann & Ellis 2000; Cole et al. 2001; Dickinson et al. 2003; 
Rudnick et alj2003Ll2006l:lFontana et : alj2004l2006l;lDrorv etalT 
12004 120051: iBorch et al.1 120061) to be estimated. Overall, if 
one assumes an st ellar initi al mass function with a power-law 
slope similar to a iSalpeteri ( 1 1955b value for masses > 1M , 
which is universa lly-applicable on galaxy-wide scales (see, e.g., 
Elmegre enl2006l for a critical discussion of this issue), the evo- 
lution of total stellar mass in the universe is reasonably well- 
reproduced by the integrated cosmic-averaged star formation 
rate (SFR) provided that one accounts for stellar mas s loss 
(ICole et al.ll20^lRudnick"eTal]l200llBorch et al.ll2006l) . 

It has also become clear that there is a bimodal distribution of 
galaxy colors at all z < 1, with a relatively narrow red sequence 
dominated by non-star- forming galaxies and a blue cloud of 
star-forming gal axies (Strate va et al.ll2001tlBlanton e t al. 2003; 
iBell et al.1 l2004t IWillmer et all l2006i see Fig. Q). More re- 
cently, is has been demonstrated that the stellar mass on the 
red sequence has built up by a substantial amount — at least 
50% of the stellar mass in presen t-day red sequence galaxies 
has come into place since z ~ 1 dChenet all 120031: belleta!.! 
2004; Fa ber et alj r2006: Bor ch et alj|2006t iBrown et alj |2006). 
Given that relatively few stars form in red-sequence galaxies, 
this build-up has been argued to be driven by the global sup- 
pression of star formation in some fraction of previously blue, 



star-forming galaxies through a variety of possible physical 
processes, e.g., galaxy merging, the suppression of star forma- 
tion in dense environments, or gas consumption by star-forming 
disks. 

Yet, this qualitative picture leaves many questions unan- 
swered. What is the contribution of dust-obscured star forma- 
tion in red galaxies to the growth of stellar mass on the red 
sequence? Can simple truncation of star formation in some 
fraction of blue galaxies provide enough stellar mass to feed 
the growth of the red sequence? The 'average' disk galaxy is 
several times less massive than typical red sequence galaxies; 
is it possible to create early-type galaxies through truncation 
of star formation in the presumably rather low-mass blue cloud 
population? 

In this paper, we present a first, crude attempt at addressing 
some of t hese questions. W e combine color-split stellar mass 
functions (IBorch et al.ll2006l) from the CO MBO-17 photomet- 
ric redshift survey JWolf etal.ll2003l 12004) with SFR functions 
derived from ultraviolet (UV) and infrared (IR) lum inosities of 
galaxies in COMBO-17 (see, e.g. jLe Floc'h et al.ll2005l for IR 
luminosity functions derived using spectroscopic redshifts and 
COMBO-17 redshifts). While both stellar masses and SFRs are 
model-dependent, depend on the assumption of a universally- 
applicable galaxy-scale stellar initial mass function (IMF), and 
suffer from considerable systematic and random uncertainties, 
such an analysis allows one to explore the growth of stellar 
mass with cosmic time and attempt to understand in which 
types of galaxies the bulk of stellar mass forms, and in which 
types of galaxies most stellar mass ends up. While necessar- 
ily qualitatitive, such an approach gives interesting insight into 
some of the basic features of the physical processes which drive 
galaxy evolution. 

The plan of this paper is as follows. In §|2] we describe 
the COMBO-17 redshift survey data, the construction of stel- 
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lar masses, and the estimation of SFR functions using Spitzer 
and COMBO- 17 data; this section includes the first estimate of 
a key observational constraint on this analysis — the estimation 
of SFR functions split by galaxy color, as a function of cosmic 
epoch. In |3] we explore the growth of the total stellar mass 
densities in red and blue galaxies, and how this growth relates to 
observed star formation in the red and blue galaxy populations. 
In $U we explore the evolution of the stellar mass functions 
of red and blue galaxies, and explore how one might expect the 
stellar mass function to evolve under various simplistic assump- 
tions. These results are discussed in ^5] In this paper, we adopt 
a Ho = 70kms _1 Mpc" 1 , 17a = 0.7 and fi m = 0.3 cosmology, and 
assume that th e distribution of stellar masses formed follows a 
Chabrier (2003) stellar IMF, averaged over galaxy scales; such 
a stellar IMF gives stellar masses and SFRs consistent with a 
iKroupal d2001l) IMF to within 10%. 

2. THE DATA 
2.1. COMBO-1 7 and stellar masses 

To date, COMBO- 17 has fully surveyed three disjoint ~ 
34' x 33' southern and equatorial fields to deep limits in 5 broad 
and 12 medium passbands. Using these deep data in conjunc- 
tion with non-evolving galaxy, star, and AGN template spectra, 
objects are classified and redshifts assigned for ~ 99% of the 
objects to a limit of mg ~ 23.5. Typica l galaxy redshift accuracy 
is fe/(l+z) ~ 0.02 (W olf etalj|2004l) . allowing construction of 
^0.1 mag accurate rest-frame colors and absolute magnitudes 
(accounting for distance and ^-correction uncertainties). Astro- 
m etric accuracy is ~ 0.1". 

iBorch etal] d2006l) estimated the stellar mass of galaxies in 
COMBO-17 using the 17-passband photometry in conjunction 
with a non-evolving template libr ary derived using the Pe- 
GASE stellar population model (see lFioc & Rocca-Volmerangd 
[19971 for a description of an earlier version of the model). The 
masses we re derived using a Kroupa et al.l d 19931) stellar IMF; 
the use of a Kroupa (2001) or 
yielded stellar masses to within < 
titatively consistent with those deriv ed using a simple color- 
stellar M/L relation dBell et al.ll2003l) . and comparison of stel- 
lar and dynamical masses for a few z ~ 1 early-type galaxies 
yielded consistent results to within their combined errors (see 
IBorch et all 120061 for more details). Random stellar mass er- 
rors are < 0.3 dex on a galaxy-by-galaxy basis, and systematic 
errors in the stellar masses (setting the overall mass scale and 
its redshift evolution) were argued to be at the 0. 1 dex level 
for galaxies without ongoing or recent major starbursts; for 
galaxies with strong bursts, masses could be overestimated by 
< 0.5 dex. In particular, the masses are derived using tem- 
plate stellar populations with an age/metallicity combination 
equivalent to roughly solar metallicity a nd ~ 6 Gyr since the 
start of star formation dBorch et alj|2006l) ; the templates do not 
evolve with redshift. While such a template set is appropri- 
ate for galaxies at z ~ 0.5-0.7 (assuming that star formation 
starts at z > 3 for most massive galaxies), galaxies at lower 
(higher) redshifts will have older (younger) underlying popu- 
lations. This leads to small, second -order offsets in stellar M/L 
at a given 'colorfl As discussed in iBell & de Jonl J200l and 
IBorch et alj d2006l see their Fig. 6), we expect that in our case 

'The first order reddening of the stellar populations as the Universe ages 
is accounted for as the aged galaxy will be fit by a redder template. This sec- 
ond order effect comes about because the older, lower redshift galaxies will 
have a somewhat larger underlying old population (carrying more mass for a 
given amount of light) than its equivalent galaxy with similar colors at high 



Chabrier! d2003h IMF would have 
10%. Such masses are quan- 



this would lead to a slight underestimate (overestimate) of stel- 
lar masses at lower (higher) redshifts at the ~ 10-20% level. 
This has the effect that the internal evolution of mass density 
within COMBO-17 will be artificially flattened. We have cho- 
sen not to apply a correction for t his oversimplificati on in what 
follows to remain consistent with Bor ch et ail ([2006). 

Stellar mass functions were derived for COMBO-17 in four 
redshift shells: 0.2 < z < 0.4, 0.4 < z < 0.6, 0.6 < z < 0.8, 
and 0.8 < z < 1. Mass-limited subsamples were constructed in 
each redshift interval: the li miting mass at eac h redshift was: 
- 1,1.8,3,6 x 1O 1O M (see IBorch et"aTI [20061 for details of 
the calculation and discussion of the limiting stellar masses). 
Mass functions (i.e., the comoving space density of galaxies 
as a function of their stellar mass) were calculated for galaxies 
above this limiting mass (i.e., no V max corrections were applied, 
as the sample is volume-limited within this redshift shell, above 
this mass limit). Schechter functions were fit to the mass func- 
tions above their respective mass limits for three different cases: 
all galaxies, blue cloud galaxies only, and red sequence galax- 
ies only. The stellar mass functions for blue cloud galaxies are 
complete well below this formal mass limit (owing to the low 
optical stellar M/L ratios of blue galaxies), therefore the extent 
to which the fitted mass function follows the data points below 
this mass limit is a test of the validity of our assumed faint-end 
slope at each redshift of interest. In each case, a different faint 
end slope was assumed: a = -1.1,-1.45,-0.7 for all, blue and 
red galaxies respectively. The se choices were m otivated by the 
z ~ stellar mass functions of Be ll et alj d2003l) : we adopt their 
stellar mass function for the local value in this work. Galaxies 
were assigned to the red sequence using the approximate cut 
U-V > 1.06-0.352z + 0.227(log 10 M*-10). The results are 
shown in Fig.Q] It is worth noting that the overall stellar mass 
functions are in excellent a greement with previous determina- 
tions (see lBorch et alj20 06), with the advantage of tracking the 
evolution of red sequence and blue cloud galaxies separately. 

2.2. MIPS 24/im data and star formation rates 

Spitzer has observed two of the COMBO-17 fields: a 1° x 
0.°5 field around the extended Chandra Deep Field South (CDFS) 
in January and February 2004 as part of the time allocated to the 
Spitzer Guaranteed Time Observers (GTOs), and a similarly- 
sized field around the Abell 901/902 galaxy cluster (A901) in 
December 2004 and June 2005 as part of Spitzer GO-3294 (PI: 
Bell). In both cases, MIPS 24/im data were taken in slow 
scan-map mode, with individual exposures of 10 s. We reduced 
the individual image frames usin g a custom data-an alysis tool 
(DAT) developed by the GTOs dGordon et alj|2005l) . The re- 
duced images were corrected for geometric distortion and com- 
bined to form full mosaics. The final mosaic has a pixel scale 
of L'25 pixel -1 and an image PSF FWHM of ~ 6". Source 
detection and photometry were p erformed using techniques de- 
scribed in lPapovich et alj d2004l) ; based on the analysis in that 
work, we estimate that our source detection is 80% complete at 
83 /iJy in the 24/im image in both the CDFS and A901 for a 
total exposure of ~ 1400 s pix" 1 . 

In order to interpret the observed 24/im emission, we must 
match the 24/im sources to galaxies for which we have red- 

redshift/earlier times. 

2 Owing to slight calibration differences between the fields, the detailed red 
sequence cut is field dependent, with the intercept at 10 10 Mq and z = being 
U- V =1.01, 1 06 an d 1.11 in the CDFS, A901 andtheSll fields respectively 
(see lWolf et aU200a for a definition of the fields). 
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FIG. 1 . — The basic ingredients of this analysis. Top panels: the evolution of galaxy colors and stellar masses over the interval 0.2 < z < 1.0, In each panel, 
the mass limit below which the stellar masses become incomplete is denoted by a thick grey line (the limit is calculated at the highest redshift in each redshift 
shell). The solid black line shows the cut used to separate red sequence and blue cloud galaxies. Central panels: The evolution of the stellar mass function. The 
stellar mass function in the red sequence (diamonds and solid line) and blue cloud (asterisks and dotted line) is shown at each redshift interval of interest. The 
z = color-split stellar mass functions are shown in each panel in gray. Uncertainties from counting statistics are shown in black; uncertainties from field-to-field 
variations are shown in gray. Lower panels: The evolution of the SFR function. Filled circles and solid lines denote the SFR function for all galaxies, asterisks 
and dotted lines blue cloud galaxies, and diamonds/dashed lines red-sequence galaxies (the blue cloud and red sequence galaxies are offset horizontally by a small 
amount for clarity). The vertical grey line shows the SFR limit, and the z = IR-derived SFR function is shown in gray in each panel for reference. Uncertainties 
are determined from the difference between our two fields with 24/im data, and our best estimate is the average of the two fields. 



shift estimates from COMBO- 17. In both the CDFS and A901 
fields, we adopt a 1" matching radius. In the areas of the CDFS 
(A901) fields where there is overlap between the COMBO- 
17 redshift data and the full-depth MIPS mosaic, there are a 
total of 3255 (4251) 24/im sources with fluxes in excess of 
83/iJy. In both fields, 70% of the 24/im sources with fluxes 
> 83/Jy are detected by COMBO-17 in at least the deep R- 
band, with mr < 26. Some 35% of the 24/im sources have 
bright niR ap < 24 and have photometric redshift z < 1 ; these 
35% of sources contain nearly half of the total 24/im flux in ob- 
jects brighter than 83/iJy. Sources fainter than niR > 24 contain 
the other half of the > 83/iJy 24/im sources; investigation of 
COMBO-17 lower confidence photometric redshifts, their op- 
tical colors, and results from other studies lends weight to the 
argument that essentially all of these sources are at z > 0.8, with 
the bulk lying at z > 1 (e.g., Le Floc'h et al. 2004; Papovich et 
al. 2004; see Le Floc'h et al. 2005 for a further discussion of the 
completeness of redshift information in the CDFS COMBO-17 



data). Given the low-confidence COMBO-17 redshifts in hand, 
we estimate at < 0.2 dex incompleteness in this 0.8 < z < 1 bin, 
and negligible incompleteness in the z < 0.8 bins. 

The goal of our analysis of the 24/im data is to obtain es- 
timates of SFR which account for both unobscured (via the 
UV) and dust-obscured star formation (via the thermal IR). Ide- 
ally, we would have a measure of the total thermal IR flux 
from 8-1000/im; instead, we have an estimate of IR lumi- 
nosity at one wavelength, 24/im, corresponding to rest-frame 
20-12/im at the redshifts of interest z = 0.2- 1. Yet, local 
IR-luminous galaxies show a surprisingly tight correlation be- 
twee n rest-frame 12-15/xm luminosity and total IR luminosity 



(e.g. ISpinodio et al.lll995klCharv & ElbazfeOOUlRoussel et all 
200U lPapovich & Belll2002l), w ith a scatter of ~ 0.15 dex. Fol- 
lowing|Papovich & Bell] ( 20021) . we choose to construct total 
IR luminosity from the o bserved-frame 24um d ata. We use 
the Sbc template from the iDevriendt et all d 19991) SED library 
to translate observed-frame 24/im flux into the 8-1000/im to- 
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tal IR luminositjQ; such a template is found to reproduce well 
the average 70 /im and 160um flux es of z ~ 0.7 galaxies with 
L IR > 10 n L Q dZheng et alj|2007al) . The IR luminosity uncer- 
tainties are primarily systematic. Firstly, there is a natural diver- 
sity of IR spectral shapes at a given galaxy IR luminosity, stel- 
lar mass, etc.; one can crudely estimate the scal e of this uncer- 
tainty by using the full range of templates fromlDevriendt et alj 
(1999), or by using templates from, e.g.. lDale et al.l (120011) in- 
stead. This uncertainty is <0.3dex. Secondly, it is possible 
that a significant fraction of 0.2 < z < 1.0 galaxies have IR 
spectral energy distributions not represented in the local Uni- 
verse: while it is impossible to quantify this error until the ad- 
vent of Herschel Space Telescope, current results suggest that 
the bulk of intermediate-high redshift ga laxies have IR spectra 
similar to galaxies in the local universe (lAppleton et aTi r2004 ; 
lElbaz et allhoollYan et al.ll2005l:IZheng et alj2007alh 



We estimate SFRs using the combined directly-observed UV 
light from young stars and t he dust-reprocessed IR emission 
of the sample ga laxies (e.g., iGordon et al.ll2000l) . Following 
iBell et all d2005l) . we estimate SFR ip using a calibration de- 
rived from Pegase assu ming a l OOMyr-old stellar population 
with constant SFR and a lChabrierld2003l) IMF: 

7/./(M yr- 1 ) = 9.8 x 10" 11 x (L IR + 2.2L UV ), (1) 
where Lir is the total IR luminosity (as estimated above) and 
Luv = ^-5iyl v ,2S00 is a rough estimate of the total integrated 
1216A-3000A UV luminosity, derived using the 2800 A rest- 
frame luminosity from COMBO- 17 Z„,2800- The factor of 1.5 in 
the 2800A-to-total UV conversion accounts for the UV spectral 
shape of a 100 Myr-old population with constant SFR, and the 
UV flux is multiplied by a factor of 2.2 before being added 
to the IR luminosity to account for light emitted longwards 
of 3000 A and shortwards of 1216A by the unobscured young 
stars. This SFR calibrati on is derived using identical assump- 
tions to lKennic utt ( 1998), and the calibration is consistent with 
his to within 30% once different IMFs are accounted for. Un- 
certainties in these SFR estimates are no less than a factor of 
two in a galaxy-by-galaxy sense, and systematic uncertainty in 
the overal l SFR scale is likely to be less than a factor of two 
(see, e.g.. lBliil200llBliir et al.ll2005l for further discussion of 
uncertainties). 

Under the assumption that the bulk of 24/im emission comes 
from star formation, not AGN activityQ we proceed to calculate 
SFR functions in a similar fashion to the stellar mass function. 
Limiting star formation rates were defined for each redshift bin: 
these are (2.6,5.2, 10.0, lS.^Moyr" 1 forz - (0.3,0.5,0.7,0.9) 
respectively. Again, the sample is volume-limited in each red- 
shift shell down to that limiting SFR. The SFR functions are 
calculated separately for all galaxies, blue cloud galaxies and 
red sequence galaxies respectively. In this contribution, we 
choose to adopt a Schechter function fit to the star formation 
rate function^, adopting a faint end slope a = -1 .45, noting that 

3 Total 8-1000/im IR lumin osities are ~ 0.3 d ex higher than the 42.5- 
122.5/im luminosities defined by Helou et al. 1 1988), with an obvious dust tem- 
perature dependence. 

4 Less than 15% of the total 24/im emission at z < 1 is in X- 
ray luminous AGN (e.g. , ISilva. Maiolino. & Gra nato 20041; IBell et~aT] [20051 
Franceschini et al. 2005; Brand et al. 2006). Although some Compton- 
thick AGN will be missed in such 



census (e.g.. | Donlev et all 2005; 
Martmez-Sansigre et al. 2006; Alonso-Herrero et al. 2006), a fraction of the 
24/im emission in galaxies with AGN is powered by star formation. Thus, 
it is likely that ~ 10% of the 24/im emission is from AGN; such a contami- 
nation is a small uncertainty compared to others inherent to this analysis, and 
does not affect our conclusions. 

5 In the local Universe, a Schechter function is a poor fit to the IR luminosity 



the total SFR density is strongly dominated by the contribution 
of blue galaxies (see Fig.[T]Q. The best-fit parameters for these 
SFR functions is given in TableQ] 

In §4, we make use of average SFRs for galaxies in given 
bins of stellar mass. These average SFRs are calculated follow- 
ing Zheng et al. (20 06), and the method ology and results are 
discussed in detail in Zheng et al. (2007b). Briefly, those galax- 
ies in a given bin in stellar mass which were not individually- 
detected at 24/im were stacked, and a total 24/im flux from 
individually-undetected sources determined. The total flux 
from individually-detected sources was then added, and the av- 
erage flux calculated. These average fluxes were converted into 
an IR luminosity, combined with UV luminosities, and were 
then used to construct SFR estimates following the above trans- 
formations. Uncertainties in average fluxes were determined by 
bootstrapping. 

3. COMPARISON OF INTEGRATED STELLAR MASS AND STAR 
FORMATION RATE DENSITIES 

It is useful to establish the basic phenomenology by com- 
paring the evolution of the integrated stellar mass density of 
the universe with the integral of the star formation rate. To 
parameterise the evolution of the total SFR, we derive the fol- 
lowing best fit to all total SFR points shown in Fig. [2] weight- 
ing equally (i.e., our measurements are included on an equal 
footing to literature measurements): log 10 psfr = — 2. 03+1.04z. 
This linear fit, shown in Fig. [2] as the solid line, appears a rea- 
sonable parameterization of the data. In order to track the frac- 
tion of star formation in red galaxies, we fit a linear function 
to the fraction of star formation in red sequence galaxies f r as 
a function of redshift: = 0.25-0.15z, defined only over the 
interval < z < 1. The blue fraction ft, = 1 — f r . Again, these 
parameterizations appear to be a reasonable description of the 
data. 

The result of this exercise is shown in Fig. [2] In the upper 
panel of Fig.|2]we show the evolution of the cosmic star forma- 
tion rate; all galaxies are shown as gray error bars (adapted from 
a compilation from Hopkins 2004) and solid circles (our mea- 
surements), asterisks denote our measurements for blue cloud 
galaxies, and diamonds denote our measurements from red se- 
quence galaxies. In the the lower panels of Fig. [2j we show the 
evolution of t he integrated stella r mass density of the Universe 
(tabulated in Borc h et aTll2006l solid circles), along with the 
integral of the total star formation rate (solid line), accounting 
for gas recy cling using the PEGASE stellar po pulation synthe- 
sis code (see lFioc & Ro cca-Volmerang d 19971 for a description 
of an earlier version of this model). O ne can clearly s ee that, 
if one assumes a universallv-applicable lChabrier! d2003l) stellar 
IMF, that the star formation rate and the growth of stellar mass 
paint a consistent picture, with around 40% of all stellar mass 
being formed at z < 1 . 

Yet, when one splits the galaxy population into red and blue 
galaxies, an interesting phenomenology emerges. It is clear that 
blue galaxies contain much of the UV+IR-inferred star forma- 

function (denoted in the lower panels of Fig.fT]by a dark gray line). Our 0.2 < 
2 < 1 measurements have a modest dynamic range in SFR, and the data can be 
well-fit with either a Schechter function or a double power-law. In this paper, 
we choose to adopt a Schechter function fit with faint end slope a = -1.45 
in order to ensure consistent behavior and extrapolations of the mass functions 
and SFR functions. Fitting with a double power-law would not affect the results 
presented in this paper to within the uncertainties. 

6 Schechter functions are fit explicitly to all galaxies and the blue galaxy 
populations; the red galaxy SFR function is estimated by subtracting the con- 
tribution of the blue galaxies from all galaxies. 
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Table 1 

Star Formation Rate Function Fits 



Zlo 


Zhi 


<f>* 


logioV'* 


a 


P 






10- 4 Mpc" 3 log 1() 


M~ l yr M yr~ l 




M Q yr~ l Mpc" 3 


All Galaxies 


0.2 


0.4 


9.8(6.0) 


1.11(27) 


-1.45 


0.020(8) 


0.4 


0.6 


10.7(1.5) 


1.20(15) 


-1.45 


0.027(4) 


0.6 


0.8 


21(5) 


1.29(12) 


-1.45 


0.065(8) 


0.8 


1.0 


17(6) 


1.42(15) 


-1.45 


0.072(7) 


Blue cloud 


0.2 


0.4 


8.3(4.9) 


1.09(34) 


-1.45 


0.016(8) 


0.4 


0.6 


8.3(2.0) 


1.16(17) 


-1.45 


0.020(6) 


0.6 


0.8 


17(5) 


1.29(15) 


-1.45 


0.055(8) 


0.8 


1.0 


15(6) 


1.45(15) 


-1.45 


0.068(7) 



Note. — Red galaxy SFR function 'fits' were determined in this paper as the overall SFR function minus the blue cloud SFR function. This was necessary owing to 
the poor number statistics of the red sequence SFR function; a direct fit frequently gave unstable results. Fit uncertainties include (rough) field-to-field uncertainties 
(half of the difference between A901 and ECDFS) and formal fitting uncertainties. Systematic uncertainties in 24/im to total IR conversion, and in the SFR calibration, 
and in the faint-end slope, are not included, and may amount to ~ 0.3 dex in the estimated SFR densities. Note that we parameterise the SFR function with a Schecter 
function. A value of Hq = 70 km s -1 Mpc -1 and a Chabrier 1 2003 ) IMF are assumed. 



tion at every epoch. If one assumes that blue galaxies will stay 
blue at all times, one can integrate the SFR to predict the growth 
of the stellar mass in blue cloud galaxies: the result is shown in 
Fig. 2 as the dotted line. It is clear that while the bulk of the stars 
are formed in blue galaxies (as traced by the UV+IR), there is 
little growth in the total stellar mass in blue galaxies with cos- 
mic time. Instead, the bulk of the growth in stellar mass is 
in red sequence (i.e., largely non-star-forming) galaxies. Blue, 
star-forming galaxies are shutting off their star formation on 
global scales and are fading onto the red sequence. 

4. WHERE IS THIS MASS TRANSFER OCCURRING? 

Further insight into the evolution of stellar mass in red and 
blue galaxies can be gained by examining the evolution of their 
stellar mass functions. While such mass functions are well- 
defined only above ~3x 10 10 M Q (the mass limit at z > 0.6), 
this exercise gives some insight into the characteristics of the 
galaxies which are transferred onto the red sequence from z = 1 
to the present-day. 

In this section, we perform a simple analysis. One starts with 
the z ~ 0.9 stellar mass function (of all, red and blue galaxies 
— whatever population is of interest). Then, one determines the 
average SFR in broad bins of stellar mass for the 0.8 < z < 1 .0 
galaxies (again, for all, red or blue galaxies), and the average 
SFR per unit stellar mass is determined as a function of stel- 
lar mass (the specific SFR dm/dt; illustrated in the left-most 
set of panels in Fig. [3). Then, the average amount of stellar 
mass growth between z = 0.9 and z = 0.7 is estimated, as a 
function of stellar mass, assuming a constant specific SFR in 
this interval and a locked-up mass fraction of 0.55 (i.e., 45% 
of initially-formed stellar mass is returned to the ISM; such 
a fraction is appropriate for a Chabrier 2003 or Kroupa 2001 
IMF for populations with age ~ 4 Gyr). Functionally, we adopt 
the simple approach of simply re-labelling the stellar mass bin 
labels of the stellar mass function with a new z = 0.7 stellar 
mass M(l + 0.55 At dm/dt). This predicted z = 0.7 mass func- 
tion (dashed line in the 2nd column of panels in Figs. [4] and |5]l 
is then compared with the observed z = 0.7 mass function (solid 
lines). The procedure is repeated between z = 0.7 and z = 0.5, 



and so on. In this way, one then can estimate (in a very crude 
fashion) the evolution of the shape of the stellar mass function 
as a function of time, and ask if this prediction resembles the 
observed evolution. 

The result for all galaxies is shown in the uppermost panel 
of Figs.|4]and|5] Given the substantial systematic uncertainties, 
and the crudeness of the analysis, the match between the pre- 
dicted evolution of the stellar mass function and the observed 
evolution is impressive. One may notice that there is a slight 
discontinuity in the behavior of the 'all' and 'red' stellar mass 
functions from the z ~ 0.3 bin to the z ~ bin. This can be seen 
also in Fig. [2] Part of this discontinuity reflects the methodol- 
ogy used for COMBO-17 mass estimation (see ^2. j} ; part is 
likely attributable to remaining cosmic variance at the low red- 
shift end of COMBO-17. 

The lower panels of Fig. |4]show the predicted evolution of 
the stellar mass function of red and blue galaxies as a func- 
tion of epoch, assuming that all mass formed in blue galax- 
ies will stay in blue galaxies for all time, and all mass formed 
in red galaxies will stay in red galaxies (i.e., no transfer of 
mass between blue and red galaxies; the assumptions used in 
constructing Fig. The low specific SFRs of red-sequence 
galaxies lead to very modest predicted mass growth, whereas a 
large amount is observed. In contrast, the higher specific SFRs 
of blue-cloud galaxies lead to the formation of a large num- 
ber of ~ 10 n M Q blue galaxies, in clear contradiction with the 
present-day blue cloud stellar mass function. 

Motivated by the overall similarity of the observed z = and 
z ~ 0.9 blue cloud stellar mass functions, we explore a sec- 
ond simple hypothesis in Fig. [5] all growth of the stellar mass 
function is added to the red sequence (i.e., the blue cloud mass 
function is not allowed to change, and the red sequence mass 
function is given all new galaxies produced in the overall mass 
functionfl This hypothesis reproduces the evolution of the blue 

7 It is also possible to carry out this exercise by adding the galaxies created 
in the blue cloud and red sequence separately. This method gives very similar 
results at the knee of the mass function, but overproduces slightly the number 
of faint red sequence galaxies. Such a discrepancy could be addressed by ar- 
guing that massive blue galaxies truncate their star formation more frequently 
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FIG. 2. — The cosmic evolution of star formation rate and stellar mass (filled 
circles and solid lines), split into contributions from the blue cloud (asterisks 
and dotted lines) and the red sequence (diamonds and dashed lines). In the up- 
permost panel, determinations of cosmic SFR are shown also in grey, adapted 
from Hopkins (2004). In the lower panels, the data points show the observed 
evolution of stellar mass as a function of redshift. The solid line shows the pre- 
dicted build-up of stellar mass, assuming the star formation history shown by 
the solid line in the upper panel (assuming for gas recycling); it is clear that the 
integral of the cosmic star formation history predicts the cosmic evolution of 
stellar mass at z, < 1 reasonably accurately. The dotted line shows the predicted 
evolution of total stellar mass assuming that all stars that form in blue galaxies 
remain in blue galaxies (i.e., if the blue cloud evolves like a closed box). The 
dashed line shows the corresponding evolution for red-sequence galaxies. In 
all panels, small horizontal offsets have been applied for clarity. 
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FIG. 3.— Specific SFR of 0.2 < z, < 1.0 galaxies. Specific SFR, defined 
as IR-derived average SFR (from stacking) per unit stellar mass, is shown for 
all galaxies, red sequence galaxies and blue cloud galaxies (top, middle and 
bottom respectively). Four different redshift bins are shown, and in each bin 
the specific SFR is shown which corresponds to roughly constant star forma- 
tion since z/ — 4 (i.e., a birthrate b of 1; grey horizontal line). Filled points 
show specific SFRs as determined for the CDFS; open points show values de- 
termined from A901. The average blue cloud galaxy has b ~ 1, red galaxies 
typically values between 1/3 and 0.1. Red sequence galaxies have significant 
star formation, as we have made no attempt to weed out star-forming systems 
with significant dust contents from the non-star-formers on the red sequence. 
The dotted line connects the average values of specific SFR. 
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FIG. 4. — Hypothesis 1. The evolution of the stellar mass function for all, red 
and blue galaxies is shown from z. ~ 0.9 to z. = (solid lines in the top, middle 
and bottom rows respectively). Dashed lines denote the predicted evolution of 
the stellar mass function, given the observed specific SFRs as given in Fig. [3] 
In this case, no transfer of mass between blue and red galaxies is permitted. 
Aitows denote the approximate limit above which mass functions are well- 
defined. 



and red stellar mass functions with striking accuracy (consider- 
ing the significant systematic uncertainties inherent to this kind 
of analysis). 

5. DISCUSSION 

In this paper, we have performed a joint analysis of galaxy 
stellar masses and star formation rates. We found an over- 
all consistency between the predicted growth of stellar mass 

than low mass galaxies. Yet, we do not overemphasize challenges at the faint 
end, as behavior of the z = 0.9 mass function at the faint end is observationally 
unconstrained. 



through star formation and the observed evolution of the stel- 
lar mass function (both in differential and integral forms). We 
showed that the observed star formation in blue cloud galax- 
ies is sufficient to dramatically overproduce the present-day 
mass budget in blue galaxies. Instead, a significant fraction of 
blue galaxies must have had their star formation suppressed by 
some physical process or processes, joining the red sequence. 
While transformation of galaxies may not be required at masses 
> 10 h Mq, large numbers of galaxies with masses < lO n M 
must be transferred from the blue cloud to the red sequence at 
z< 1. 



Star formation and the growth of stellar mass 



7 



-2 

-3 

I -1 

o 

„ -2 
u 

°- -"5 

^ -4 

2 -1 

a* 

o _ 2 
-3 
-4 



z ~ 0.9 


z ~ 0.7 


z ~ 0.5 


z ~ 0.3 


z - 0.0 






^\ 
















^Red ^\ 
















\ 

V 


\\ 
\\ 



10 11 12 10 11 12 10 11 12 10 11 12 10 11 12 
log, M./M 



FIG. 5. — Hypothesis 2. This figure is formatted similarly to Fig. [4] except 
that in this case all newly-formed galaxies are assigned to the red sequence. 
That is, the blue cloud is not permitted to grow, and the global growth in galaxy 
numbers at a given mass are assigned to the red sequence. 



5.1. The empirical basis for this paper 

While the analysis presented in this paper is clearly model- 
dependent, the main results are qualitatively robust. The em- 
pirical basis for this statement is clear. The rest-frame op- 
tical luminosity density in blue cloud galaxies decreases by 
~ 1 mag since z = 1 ; the change is to fir st order a fading 
in L* without accompanying chan ge in (jf dLillv et alJll995b 
IWolf et al.ll2003tlWillmer et alj|2006l) . In this interval, the rest- 
frame optical luminosity density of red sequence galaxies stays 
roughly constant; the sense of this change is that L* fades by 
~ 1 mag since z = 1 , accompanie d by an increasing <jf to keep 
the luminosity density constant ( Chen et all 120031: iBell et all 
l2004t iFaber et alj 120061: iBrown et alj 120061) . The colors of 
both red and blue galaxies beco me redder with increasing time 
dBell et al.ll2004t lBlantonll2006l) : it is impossible to avoid some 
increase in mass-to-light ratio with such aging dBell & de Jongl 
l200lMKauffmann et al.ll2003HlBell et alj|2003l) . Thus, the stel- 
lar mass function of blue cloud galaxies must be approximately 
unchanging, while the stellar mass function of red sequence 
galaxies must undergo sign i ficant growth, at least around L * 
(see, e.g.. lBundvetal1 l2"006; Bla ntonll2006l:lBorch et al.ll2006h . 
On the other hand, UV, line emission and thermal-IR surveys 
all demonstrate convincingly that the bulk of star formation is 
in blue cloud, primarily disk-dominated gala xies ( Flores et alj 
1999 :IZheng et al.ll2004 iHammer et alJl2005l: iBeil et alJl2005t 
Melbour ne et alj 120051) . The IR emission, in particular, ar- 



gues for significant extra dust-obscured star format ion above 
and beyond what is seen in the local Univers e (Elbaz etalj 
1999 :ICharv & Elbaz|200"ll;lTakeuchi. Buat & Burgarellal2005t 
Le Floc'h et alj|2005l) . suggesting typical star formation rates 
for massive disks of > 10M Q yr (i.e., doubling times of or- 
der a Hubble time; Hammer et al. 2005). Thus, the massive 
blue galaxies contain the bulk of the star formation; since one 
is not allowed to change significantly the number of massive 
disk galaxies (because of the luminosity function and color con- 
straints), these rapidly-growing blue sequence galaxies must 
be balanced by a loss to the (also rapidly building-up) red se- 
quence. In this way, despite the many model dependencies 
— the choice of bolometric correction for 24/im-derived star 
formation rates, calibrations of star formation rates and stellar 



mass-to-light ratios, choices about faint-end slopes, choice of 
IMF, etc. — it is clear that the overall picture presented in this 
paper is rather robust. 

5.2. 3 x 10 10 Mq: More than just the crossing point of mass 
functions ? 

A particularly interesting result of this analysis is that the 
mass scale at which the bulk of the red sequence growth takes 
place is approximately the same as the mass scale at which the 
bulk of star formation occurs on the blue cloud, at between 
1O 1O M and 10 n M o ; i.e., at - 3 x 10 10 M Q . It was not ob- 
vi ous at all t hat th is needed to be so. We ourselves argued 
in IBell et alJ d200l that there were an insufficient number of 
massive disk galaxies with which to feed the red mass function 
growth: faced with this mismatch, we appealed to widespread 
and frequent merging of low mass ~ 1O 1O M galaxies to pro- 
duce more massive red sequence counterparts. Such a view 
is manifestly mistaken, however. Given the results from ISO 
and Spitzer, it is now clear that massive disk galaxies host a 
large amount of star for mation, sufficient to double their mass 
in a Hubble time or les s dHammer et al.ll2005llMelbourne et alj 
20051 IBell et~aT]|2005h . Such intense star formation, obscured 
by dust, is sufficient to dramatically overproduce the number 
of local massive disks. This shifts our perception of the blue 
cloud significantly. A large fraction of massive blue galaxies 
are required to quench their star formation and join the red se- 
quence; this loss of massive blue galaxies is compensated for 
by growth of lower-mass blue cloud galaxies through intense 
star formation. 

In the past, much significance has been attached to the tran- 
sition in galaxy properties see n at ~ 3 x IQ 10 M P ,, prese nted 
first (in a compelling form) by Kauffman n et alj ( 120 03b). In 
its original incarnation, this transition is defined as the mass 
at which blue galaxies and red galaxies are equally common, 
i.e., where the blue and red galaxy stellar mass functions cross. 
Subsequent work has found other galaxy properties to change at 
this transition mass, e.g., AG N activity is much more common 
above this mass than below it (Kauffmann et al. 2003c). Recent 
works have probed the evolution of the cross-over of the red and 
blue mass functions out to z = 1, finding an order of magnitude 
decrease in the cro s s-over point from z = 1 to the present day 
dBundy et al.l l2006t iBorch et alJ 12006). We take the position 
that this cross-over point is to a certain extent a historical hap- 
penstance; indeed, it is observe d to be a strong fu nction of envi- 
ronment in the local Universe dBaldrv et alJ l2006). Instead, we 
attach much more significance to the cut-off of the blue cloud 
stellar mass function, starting at ~ 3 x 10 10 Mq. This cut-off, 
which on the basis of this analysis we interpret as an important 
scale for global suppression of star formation and transforma- 
tion to red sequence galaxies, is not a strong function of redshift 
or environment. 

5.3. Expectations of this picture 

There are three predictions/expectations of this approach 
which can be compared with independent analyses. Firstly, it 
is increasingly apparent from previous analyses as well as our 
own that massive red sequence galaxies with M* > 10 n M Q 
should have older stellar p opulation than those at 3 x 10 10 



should have older stellar population than those at j x 1U 
(see. e.g..Uuneau et alJfcOQSl: ICimatti. Daddi. & Renzinill2006l; 



iPanter etaT. 20061) . On the basis of the observed b <~ 1 of 
blue cloud galaxies (galaxies whose SFR has not yet been sup- 
pressed) and the observed growth of the red sequence, one 



8 



Bell et al. 



might expect that the average age of a red sequence galaxy 
with 3 x 10 10 M Q should be < 8 Gyr. At masses > 2 x 10 U M Q , 
little growth is observed; one would expect m ost star forma- 
tion to happen at z > 1 (as is indeed observed; Pap ovich et alJ 
200d iDaddi et al.ll2005h . giving typical 'ages' of ~ 8- 12 Gyr. 
This compares favorably with observati ons of the stellar pop - 
ulation ages of ea rly-type galaxies fr om Thomas et alJ J2005), 
discussed also by [Thomas & Daviesl d2006l): galaxies with ~ 
3_xl0 10 Mo have ages between 5 and 9 Gyr dThomas & Daviesl 
2006, adopting the age scatter of ~ 0.25f, where t is the age, 
from Thomas et al. 2005), w hile galaxies with masse s > 2 x 
1O U M have ages 9- 13Gyr dThomas & Daviesll2006l) . 

Secondly, on a related theme, one expects also a relatively 
prominent population of galaxies undergoing truncation of their 
star formation at masses between 10 10 Mq and 10 n M Q . Such 
a population will be characterized by a rapidly-declining or 
truncating star formation rate; depending on the rapidity of 
the truncation, and on the distribution of dust in the galaxy, 
one is likely to see some spectral evidence of the rapid de- 
crease star formation rate (e.g., enhanced Balmer lines and 
weak/non existent emission lines) . In this context, the obser- 
vations of Le Borg ne et al.l d2006l) are particularly interesting. 
They study galaxies with M > 1.7 x 1O 1O M at 0.6 < z < 1.2 
and the present day using the Gemini Deep Deep Survey and the 
Sloan Digital Sky Survey, finding a decrease in the "H<5 strong" 
fraction from ^50% at z ~ 1 .2 to a few percent at present. Ex- 
tension of this methodology to larger samples (allowing split- 
ting into different mass bins) will allow a more detailed investi- 
gation of the redshift evolution, environmetal dependence, and 
mass dependence of the truncation of star formation in galaxies. 

Finally, depending on transformation mechanism, one may 
expect quite a bit of environmental dependence in galaxy ages; 
exploration of this issue is well beyond the scope of this work. 
Qualitatively, however, given that in this picture a large frac- 
tion of current red galaxies are generated through suppression 
of star formation in galaxies which were blue at z ~ 1, and that 
red galaxies currently live in reasonably dense environments, 
one expects that a large number of massive blue z ~ 1 galax- 
ies will be relatively strongly clustered — clustered much more 
strongly than present-day massive blue galaxies. Such an ex- 
pectation is spectacularly borne out by the data: analyses of 
both the DEEP2 and VVDS galaxy evolution surveys demon- 
strate an abun dance of unusually strongly-clustered m assive 
blue galaxies dCooper et alJ 120061 ; ICucciati et a n i2006h . For 
workers wishing to understand the driving forces of early-type 
galaxy evolution, these galaxies merit special attention: anal- 
yses of their environment, morphological and structural prop- 
erties, stellar populations, and supermassive black hole content 
will shed light on these 'proto-early-type' galaxies. 

6. CONCLUSIONS 

In this paper, we have presented a joint analysis of galaxy 
stellar masses and star formation rates from the COMBO- 17 
photometric redshift survey combined with Spitzer 24/.tm data. 
Stellar masses were estimated using 17-passband optical spec- 
tral energy distribution fitting, and star formation rates were de- 
termined from the combined ultraviolet and IR fluxes of galax- 
ies. We estimated SFR functions for all galaxies in 4 redshift 
bins, and presented for the first time SFR functions split into 
contributions from blue cloud and red sequence galaxies. 

When standard stellar mass and SFR calibrations were used, 
we found an overall consistency between the observed growth 
of stellar mass since z = 1 and the integral of the cosmic SFR. 



Yet, the bulk of SFR occurs in blue cloud galaxies; if all stars 
either already in or formed in blue galaxies since z = 1 were to 
end up in present-day blue cloud galaxies, the total stellar mass 
budget in blue cloud galaxies today would be dramatically over- 
produced by a factor of two. Instead, the stellar mass density in 
red sequence galaxies grows steadily from z = 1 to the present 
by approximately the amount required to balance the 'overpro- 
duction' of stars in blue cloud galaxies. Thus, a large fraction 
of blue galaxies must have their star formation suppressed by 
some physical process or processes, joining the red sequence. 

We explored the specific SFRs (i.e., SFR per unit stellar 
mass) as a function of galaxy mass. These specific SFRs were 
calculated using stacking of 24/im data on the positions of the 
galaxies of interest, in combination with their rest-frame UV 
luminosities from COMBO-17. Using these specific SFRs, we 
'predicted' the evolution of the stellar mass function from one 
redshift bin to the next. We found that the evolution of the 
overall stellar mass function is reasonably well reproduced us- 
ing this approach. Yet, if one assumes that all stars already in 
or formed in blue galaxies since z < 1 stay in blue galaxies, one 
dramatically overproduces the number of present-day massive 
blue galaxies and underproduces the growth of red sequence 
galaxies at < 10 u Mq. This predicted excess of massive blue 
galaxies is approximately sufficient to feed the growth of the 
red sequence at < 10 u Mq. Thus, not only must a significant 
fraction of blue star forming galaxies suppress their star forma- 
tion between z = 1 and the present day, but also the observations 
require that the suppression mechanism or mechanisms be par- 
ticularly effective at ~ 3 x 10 10 M Q . 
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Spitzer Space Telescope Fellowship Program, through a con- 
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